The question of whether DNA is intrinsically conducting or 
not is still a challenge. The ongoing debate on DNA molecules 
as an electronic material has so far underestimated a key dis- 
tinction of the system: the role of base pairing in opposition 
to correlations along each chain. We show that a disordered 
base paired double-chain present truly or, at least, effectively 
delocalized states. This effect is irrespective to the sequencing 
along each chain. PAC(s) numbers: 72. 15.Rn,72.80.Le, 73.22.- 
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The convergence of two scientific branches - design of 
nanometric low dimensional systems in condensed matter 
and controlled growth of nucleotide sequences in molec- 
ular biology - leads to the quest of DNA as a new na- 
noelectronic material [1]. The main question concerning 
the intrinsic conductivity remains unsolved. Experiments 
on DNA conductivity are very controversial: metallic [2] , 
semiconducting [3] , insulating [4] and even superconduct- 
ing [5] behaviors have been reported. A recent work, al- 
though, has shown that electrons or holes are responsable 
for electrical current in DNA [2] and a biological conse- 
quence is the fact that a mechanism for sensing dam- 
aged bases may explore the long range electron migra- 
tion along the molecules [6] . These experiments are very 
complex due to the influence of the local environment 
such as counterions, contact resistance, thermal vibra- 
tions and even sequence variability, which are difficult to 
control in non designed samples [7]. 

Here we focus on an intrinsinc factor: the role of base 
pairing in opposition to correlations along each chain as a 
fundamental mechanism governing the conductivity. The 
study of correlations effects on disorder has already a long 
history. In summary, disordered one dimensional sys- 
tems should show only stricly localized eletronic states 
[8] in absence of electronic interactions. However, when 
a short-range correlation is imposed on disorder, delo- 
calized states arise, like in a random dimmer correlated 
one-dimentional chain [9]. On the other hand, DNA can 
be enginnered to almost any imaginable sequence, there- 
fore a completely random sequence is an important limit- 
ing case. Some authors claim that long range correlations 
could be present in some genes [10,11] but a further claim 
that such correlations could lead to delocalization need 
further study [12]. 

The theoretical studies of the electronic properties of 
DNA range from strictly one-dimensional tight-binding 
chains [13] [14], up to involved ah initio and density- 
functional methods [7] [15]. Although partially success- 
ful, there are severe limitations in both approaches. One- 
dimensional chain models deal with effective sites instead 
of a double-chain structure. Therefore, the base pairing 
is not properly taken into account. On the other hand, 
involved density functional calculations, although giving 
usciful insights on environment influence [15], have to be 
limited to a reduced number of model DNA molecules. 
Nevertheless, numerous models in the literature are in- 
termediate to both limits, like the charge ladder model 
[16] in which two ordered tight-binding chains are consid- 
ered together with a Coulomb repulsion between bases. 



Another approach catching the double chain character of 
the DNA molecule focus on short ordered molecules with 
a backbone chain [17], hindering a clear study of the base 
pairing effect. 

In the present work, we show that DNA-like base pair- 
ing leads to delocalization in double chains. Two funda- 
mental aspects are discussed. First, there is indeed a true 
localization-delocalization transition for certain parame- 
ters ranges. Secondly, we show that at least a very effec- 
tive delocalization, in the microelectronic length scale, 
is induced by base pairing for a wide parameter range 
that is compatible with the DNA electronic structure. 
Furthermore, it is important to notice that this delocal- 
ization is irrespective of correlations along the DNA se- 
quencing. We use a double-chain nearest-neighbor (intra 
and inter-chain) hopping tight-binding model to describe 
the system. In this approximation, the Hamiltonian can 
be writen as: 

N/2 2 

H = Y,i^[^ij\i,3)(jA+V\i + l,3)(j,i\ + V\i-l,3)(j,i\ 

+V'\i,j + l)(j,z|5,,i + V'\i,3 - l)(j,zl<5,,2] (1) 

where £i,j is the (i,j) site energy, V is the intra-chain 
hopping parameter, V' is the inter-chain hopping param- 
eter and N /2 is the number of base pairs (N is the total 
number of sites). 

The double chains are constituted by four different 
sites, representing the nucleotides A, C, G and T. These 
sites are randomly assigned in the first chain with equal 
concentrations on average, while the sites of the second 
chain follow the DNA pairing. It should be noticed that, 
unlike the costumary single-chain correlation (intra-chain 
correlation) [13] [14], this is a base pairing correlation, 
which should necessarily be taken into account in or- 
der to catch the DNA key features. The inset in Fig. 
1(a) shows a particular double-chain configuration. The 
present model is completely general but will be tied to 
the quest of the electronic properties of DNA by means 
of a proper parametrization. We use the following site 
energies: £a = 8.24 eV, Et = 9.14 eV, ec = 8.87 cV and 
EG = 7.75 eV, according to values suggested in the litera- 
ture [18]. The hopping parameter for intra chain adjacent 
nucleotides is taken to be = 1.0 eV. The inter chain 
coupling, V , is also treated as a hopping integral. Ini- 
tially, as a model of DNA molecules we consider V < V, 
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as expected [4,19]. All important consequences of a wider 
parameter range will be discussed at the end. To avoid 
spurious effects due to a particular configuration, an av- 
erage over 20 double-chains is always undertaken. 

The degree of localization of a state is given by the 
participation ratio (PR) [20], defined, in tight-binding 
approximation, by: 

where aij is the wave function amplitude in the site. 
The PR is close to zero for localized states for ^ oo 
and reaches the maximum value of 2/3 for delocalized 
states in one dimensional ordered systems [21]. 

In order to reveal the effect of base pairing, we compare 
base paired double chains with the uncorrelated counter- 
parts. Fig. 1(a) shows the density of states (DOS) for 
coupled double chains with 1250 base pairs. The hopping 
parameters here are V = 1.0 eV [18] and V = 0.5 eV. 
There is a no qualitative difference between base paired 
(indicated by an arrow) and uncorrelated double chains: 
both show the expected van Hove singularities for a dou- 
ble chain, which are not destroyed by the alloy disorder, 
although the peaks are more pronounced for the base 
pairing case. On the other hand, the base pairing intro- 
duces a dramatic change in the degree of localization of 
the states, as can be seen in Fig. 1(b). The base paired 
chains show two bands of effectively delocalized states 
(indicated by arrows) [22]. In the limit of V = eV 
the system reduces to two completely random chains. As 
soon as an inter chain coupling is turned on, together 
with the base pairing, the delocalization mechanism is 
induced (as shown by the raise of the PR values). It 
would be interesting to verify the limit in which the in- 
ter chain coupling dominates over the intra chain one, i. 
e., V > V. It is worth mentioning that bands of truly 
delocalized states appear due to base pairing if the inter 
chain hopping is greater than the intra chain one. The 
inset of Fig. 1(b) shows the participion ratio for a base 
paired and a uncorrelated double chain with N/2 = 2000 
base pairs with the following sites energies ea = 10.5 eV, 
£T = 10.3 eV, sc = 9.5 eV and ec = 9.7 eV. These 
values are close to the DNA parameters [18] but chosen 
to be symmetric relative to a reference level at 10.0 eV. 
The hopping parameters are: V = 2.0 eV and V = 3.0 
eV. The base paired case shows two bands with the PRs 
given by smooth curves reaching 2/3, the expected value 
for completely delocalized states in ordered chains [21]. 
The two peaks with even higher PRs are fingerprints of 
edge states coiiplcd to the chain continuum. 

The degree of localization of an electronic state can 
be visualized by the extension of the corresponding wave 
functions. We choose a single disorder configuration of 
the system discussed in Fig.l and look at eigenstates 
with energies near the PR maximum around 10.7 eV. 



Fig. 2(a) depicts a wave function spreadcd out along a 
entire chain of the base paired system; while the inset 
shows a strongly localized state at the same energy for 
the corresponding uncorrelated double-chain. 

Up to now we have focused on the extreme limit of 
a completely random DNA sequencing. For a random 
Poly(A)-Poly(T) double strand, the delocalization mech- 
anism is even more important. In Fig. 3 we show the 
PR for base paired (indicated by an arrow) and unpaired 
Poly(A)-Poly(T) double chain (base paired) N/2 = 3000 
base pairs long. The wave function of a typical delocal- 
ized state of the base paired system is shown in the inset 
[23] . Effective delocalization has been verified up to 5000 
base pairs long double chains (not shown here). 

A crucial test for the delocalization is the evolution of 
PR with the length of the system. Fig. 4 shows a map of 
the average maximum PR as a function of double chain 
length - up to 3000 base pairs - and inter chain hopping 
for the completely random sequencing limit. The intra 
chain hopping is fixed at V = 1.0 eV. Fig. 4(a) illus- 
trates base paired cases, while Fig. 4(b) is a mapping 
of the double chains in absence of base pairing. The un- 
coupled chains limit, V = 0.0 eV, is stricly identical in 
both systems, as expected. As soon as the coupling is 
turned on, the base pairing leads to higher PR values for 
any value of inter chain hopping. In order to better il- 
lustrate this behaviour we mark one of the contour levels 
as a guide for the eyes. The horizontal lines represent 
the DNA-like case discussed in Figs. 1 and 2. The con- 
tour level crosses the horizontal line for the base paired 
case. Fig. 4(a), at a chain length 3 times longer than 
for the uncorrelated case. Fig. 4(b). For stronger inter 
chain couplings, this length ratio - between base paired 
and uncorrelated chains - becomes continuously higher. 
This suggests a crossover from localized to truly delocal- 
ized systems, as pointed out for the case shown in the 
inset of Fig. 1(b). Such localization/dclocalization tran- 
sition can be verified quantitatively by investigating the 
participation number (PN) as a function of chain length. 
[24] 

The PN, given by multiplying the PR by the chain 
length, should increase with the chain length as PN = 
N X PR oc N'^^ [24]. In summary, for a localized state 
D2 = 0, i.e., PN is independent in respect to N. For 
delocalized states D2 = 1, in the present case, while 
< < 1 occurs for wave functions that are effectively 
delocalized, but occupy only a fraction of the system in 
the thermodynamic limit. In Fig. 5 the average maxi- 
mum s for various double chains are depicted as a func- 
tion of ln(N/2). The dashed curve corresponds to base 
paired double chain (DNA-like) discussed in Figs. 1 and 
2. The continuous line coresponds to the uncorrelated 
case (also shown in Figs. 1 and 2). For the lengths inves- 
tigated, we clearly see the localized behavior in the un- 
correlated case ( approximately constant PN). The base 
paired case (dashed curve) shows an effectively delocal- 
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ized character. It is interesting to notice that truly delo- 
cahzed case, shown in the inset (see also the inset in Fig. 
1(b)): the exponential growth in the present represen- 
tation corresponds to the linear behavior for delocalized 
states. 

In conclusion, we investigate the effect of base pairing 

on the transport properties of DNA-likc systems. Delo- 
calization in coupled one-dimensional chains has been re- 
ported only for a very particular case [25]: an odd number 
of uncorrelated chains and at a single critical energy. In 
the present work we now impose a base pairing that mim- 
ics the key feature of DNA molecules and entire bands of 
delocalized states arise. This delocalization is effectively 
originated by the base pairing, since the system has an 
even (2) and not an odd number of chains [25]. On the 
other hand, a localization-delocalization transition has 
been discussed for double chains, but with intra- chain 
correlated disorder [26] but in absence of base pairing. 
Both results do not apply to a DNA like struture. In our 
model the base pairing induces dramatic effects on the 
degree of localization: base pairing leads to an effective 
localization-delocalization transition irrespective to the 
sequencing. A recent report presents experimental evi- 
dences of delocalization of injected charges in DNA [27] 
that could be related to base pairs coupling. The present 
results suggest that DNA is intrinsically a promising elec- 
tronic material and the hindrance to DNA nanolectronics 
is solely of technological nature. Finaly, we show strong 
evidences that base pairing also leads to a bona fide lo- 
calization delocalization transition at certain parameter 
ranges. 
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(a) Density of States (DOS) of a disordered double 
chain with (dashed) and without (continuous line) base 
pairing. The length of the chain is 1250 base pairs. In- 
set: schematic representation of a double chain with C-G 
and A-T base pair inter-chain correlation (b) PR for the 
double-chains in (a): with (indicated by an arrow) and 
without base pairing. Inset: The same but for different 
parameters (see text) 

Wave function for an cigenstate near the PR maxi- 
mum in Fig. 1(b) (Double-chain with base pairing) Inset: 
a localized state for the uncorrelated case at the same 
energy. 

PR for a Poly(A)-Poly(T) double chain N/2 = 3000 
base pairs long. Base paired (indicated by an arrow) and 
uncorrelated cases. Inset: a delocalized wave function for 
the base paired chain 

Conture plots of the average PR maximum as a func- 
tion of the chain length and inter chain coupling, (a) base 
paired double chains, and (b) double chains without base 
pairing. For details see text. 

Participation number function of double chain 

length: with (dashed curve) and without (continuous 
curve) base pairing. Same parameters as in Fig. 1. Inset: 
the same plot for the system in the inset of Fig. 1(b). 
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